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We consider the role of static disorder in the spin sector 
of the one- and two-channel Kondo models. The distribution 
functions of the disorder-induced effective energy splitting be- 
tween the two levels of the Kondo impurity are derived to 
the lowest order in the concentration of static scatterers. It 
is demonstrated that the distribution functions are strongly 
asymmetric, with the typical splitting being parametrically 
smaller than the average rms value. We employ the derived 
distribution function of splittings to study the temperature 
dependence of the low-temperature conductance of a sample 
containing an ensemble of two-channel Kondo impurities. The 
results are used to analyze the consistency of the two-channel 
Kondo interpretation of the zero-bias anomalies observed in 
Cu/(Si:N)/Cu nanoconstrictions. 
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I. INTRODUCTION 

The Kondo effect, that is the low-temperature screen- 
ing of dynamical quantum defects in metals by band 
electrons, has been extensively studied during the past 
thirty years and is by now well understood (see, e.g., 
Refs. [0,0). The continuing interest in the problem is 
motivated chiefly by the search for novel realizations of 
the effect p] . For example, local non- Fermi-liquid ground 
states have been predicted theoretically Q for certain 
types of dynamical defects coupled to band electrons. 
The behavior of these impurities is sometimes invoked 
in efforts to understand the non-Fermi-liquid behavior of 
strongly correlated systems, such as heavy fermion ma- 
terials and high-temperature superconductors H. 

One of the main difficulties encountered in the inter- 
pretation of experimental data from novel Kondo systems 
is the fact that defects with internal degrees of freedom 
very seldom represent the only type of disorder in the 
system. More often, a considerable number of random 
static defects affecting band electrons are also present. 
Scattering of electrons on this static disorder may alter 
the experimental signature of the Kondo effect [p|-fi2[, 
sometimes masking genuine non-Fermi-liquid behavior, 
or, vice versa, possibly mimicking it in Fermi-liquid sys- 
tems |3). 

The Kondo effect requires that at least two internal 
states of the impurity be degenerate or, at least, that 
their energy difference be much smaller than the Kondo 
temperature Tk. Barring the cases of accidental degen- 
eracy between the states of the dynamical impurity, de- 



generacy occurs as a consequence of a symmetry, e. g., 
invariance under time-reversal transformation in the case 
of the magnetic one-channel Kondo effect. Accordingly, 
the possible types of disorder can be separated into two 
classes, depending on whether disorder destroys the rel- 
evant symmetry of the Hamiltonian. 

A typical example of the first class is a dilute magnetic 
alloy with a finite concentration of non-magnetic defects. 
Assuming that the charge state of the magnetic impu- 
rity does not change as a result of the interaction with 
conduction electrons (thus excluding the mixed-valence 
regime of the Anderson model), potential scattering on 
static disorder does not involve the degree of freedom of 
conduction electrons which is coupled to the magnetic 
impurity - their spin. The Hamiltonian remains invari- 
ant under time reversal, and the spin states of the impu- 
rity are degenerate even in the presence of disorder. The 
Kondo temperature Tr-, on the other hand, is affected 
by potential scattering of electrons. The nature of the 
ground state in systems of this type has been studied in 
Refs. 110,111. 

Another example of such symmetry-preserving disor- 
der is given by spin-orbit scattering lljJTS] in magnetic 
Kondo systems. The corresponding Hamiltonian is also 
invariant under time reversal, and Kramers' theorem en- 
sures that each orbital state of conduction electrons is 
doubly degenerate, so their coupling to magnetic impu- 
rities does not lift the degeneracy of impurity states. 

An entirely different situation is encountered when 
scattering on static defects breaks the relevant symme- 
try. Kondo coupling between the band electrons and the 
dynamic defect then leads to symmetry-breaking contri- 
butions to the self-energy of the dynamic defect. The 
frequency-independent part of these contributions can 
be reinterpreted as an extra term in the bare Hamilto- 
nian of the dynamic defect, and the difference between its 
eigenvalues (the energy difference between the "up" and 
"down" states of the defect) is an effective splitting A. 
Being induced by random scattering, the splitting itself 
is a random quantity. This type of model was studied, for 
example, in Ref . iQ , where the average splitting induced 
by broken time-reversal invariance due to the combined 
effect of random spin-orbit scattering and weak magnetic 
field was computed. A similar model has been encoun- 
tered in the study of internal magnetic field distributions 
in spin glasses Q . 

The present study of models of this type has been 
motivated in part by the discussion in Ref. |ll| of dif- 
ferent theoretical interpretations of zero-bias anomalies 



in Cu/Si:N/Cu nanoconstrictions [0|. The zero-bias 
anomalies first reported in Ref. |17|] were observed in 
nanoconstrictions formed by etching a bowl-shaped cav- 
ity in an insulating Si3N4 substrate before covering both 
sides with vacuum-deposited Cu films. The anomalies 
are characterized by a \fV dip in conductance G at small 
bias V ^ a corresponding vT temperature dependence of 
conductance at 1^ = 0, and, more generally, a scaling 
function of the form G (V, T) - G (0, T) = T^/^P {V jT) 
where F (a; ^ 1) a x^l"^ . These features were interpreted 
in Refs. ||l^,|l8| (see also Ref. [ll^ for additional exper- 
imental results and Ref. pOf for a comprehensive re- 
view) as consistent with the scaling properties at low 
temperatures of the two-channel Kondo model [Q . The 
observed absence of Zeeman splitting led to a conclu- 
sion p7[ that a non-magnetic realization of the two- 
channel Kondo model of the type suggested by Vladar 
and Zawadowski [El| might be responsible for the ob- 
served anomaly. 

The two-channel Kondo model, first proposed in 
Ref. Q (see also Ref. [^ for an extensive review) to clas- 
sify magnetic properties of rare earth materials, is char- 
acterized by doubling of the degrees of freedom of con- 
duction electrons as compared to the usual one-channel 
case, while the dynamic impurity is still a "spin-up, spin- 
down" doublet. In other words, each orbital state of 
conduction electrons acquires, in addition to its spin, an 
extra label, "flavor", which is silent in the sense that 
the scattering on the dynamic impurity conserves the 
flavor quantum number. Even so, the strongly corre- 
lated ground state of this model has been predicted to 
exhibit unusual and rather distinctive scaling properties, 
markedly different from the Fermi-liquid-like ground- 
state properties of ordinary Kondo impurities. 

In the original model proposed in Ref. 0], the flavor 
degrees of freedom were constructed out of different an- 
gular momentum states of conduction electrons. Subse- 
quently, it was proposed in a series of papers by Vladar 
and Zawadowski [gl|] that an effective two-channel Kondo 
model may emerge in an entirely different context, where 
the roles of orbital angular momentum and spin of con- 
duction electrons are interchanged. The role of a dy- 
namic impurity in such non-magnetic realizations of the 
Kondo effect is assumed to be played by a two-level sys- 
tem (TLS) - an atom or a group of atoms tunneling be- 
tween two nearly degenerate states. If transitions be- 
tween the two states of the TLS involve a transfer of 
charge, the transition amplitude becomes dependent on 
the density of conduction electrons via the Coulomb in- 
teraction 1^ . The parity of electronic states with respect 
to the center of the spatially extended defect becomes the 
active degree of freedom - "pseudospin" . The physical 
spin assumes the role of the silent "flavor" degree of free- 
dom, providing two independent channels (in the absence 
of spin scattering) for the screening of the dynamic im- 
purity y,|2^. Such TLS may be formed accidentally in 
a strained glassy material Q , or as a result of a Jahn- 
Teller effect [M, and TLS have also been conjectured to 



occur at interfaces |£5|. Very recently, the non- Fermi- 
liquid properties of the ground state in this model have 
been invoked in the study of dephasing rate of conduction 
electrons in disordered metals P6|] . 

The degeneracy between the two states of the TLS 
based on pseudospin symmetry is a crucial precondition 
for two-channel Kondo screening and the formation of 
the non-Fermi-liquid ground state. In practice this de- 
generacy is almost always expected to be lifted because 
the pseudospin, corresponding to parity about the cen- 
ter of the TLS, is not in general a conserved quantity, 
and, in particular, is not conserved by ordinary potential 
scattering. This feature of the orbital two-channel Kondo 
model has to be contrasted with magnetic Kondo models 
where the relevant symmetry is time-reversal invariance, 
which is broken only in special circumstances, e.g. by an 
applied magnetic field or magnetic disorder. Therefore, 
an analysis of the magnitude of disorder-induced split- 
tings of two-level systems is essential in evaluating the 
consistency of the two-channel Kondo interpretation of 
the zero-bias anomalies of Refs. [ pj|p^ . 

Most of the previous theoretical work on this sub- 
ject [^ 14 1 has been concentrated on computing the sec- 
ond moment of the random splittings. A calculation of 
<^A^) induced by white-noise potential scattering was re- 
ported in Ref. ||ll[. It was argued there that even small 
amounts of disorder may lead to large splittings between 
the energy levels of the TLS, thus effectively stopping 
the Kondo screening at temperatures higher than Tr-. 
However, the distributions of splittings tend to be very 
asymmetrical so that their moments are not representa- 
tive of the typical values. Moreover, the knowledge of 
the full distribution function is necessary to understand 
how the splittings of an ensemble of defects may affect 
the scaling behavior of conductance. 

It should be remarked that there is no direct evidence 
for the existence of TLS in the nanoconstrictions studied 
in the experiments of Refs. p? 19|. It is precisely the 
match between the experimental scaling of conductance 
and that predicted theoretically for two-channel Kondo 
systems that is the main argument in favor of the two- 
channel Kondo interpretation of the data. The theoreti- 
cal scaling functions used for this purpose in Refs. [^ 18 
were derived under the assumption that no disorder other 
than the TLS themselves is present, and thus it is of con- 
siderable interest to understand how these scaling func- 
tions may be changed by realistic amounts of static dis- 
order. 

In this paper we study distribution functions of split- 
tings for two models: (i) isotropic magnetic Kondo impu- 
rities in a spin glass environment, and (ii) atomic TLS in 
an environment of static defects inducing potential scat- 
tering. In both cases the disorder is modeled by an ar- 
ray of randomly located point scatterers. The magnetic 
model may be realized if a dilute solution of weakly cou- 
pled magnetic impurities undergoes Kondo screening in a 
spin-glass environment, formed, for example, by a more 
concentrated solution of more strongly coupled magnetic 



impurities. The result for the distribution of splittings 
in this case reproduces the distribution of internal fields 
in spin glasses derived earlier in Ref. [|l6[. We also ana- 
lyze the effects of higher-order terms in the Kondo cou- 
pling which cannot be reduced to RKKY-type expres- 
sions. These terms are shown to lead to a finite renormal- 
ization of the small-A portion of the distribution func- 
tions, while their effect on the large-A tail is negligible. 
A similar analysis is performed for the distribution of 
splittings of atomic TLS. The magnitude of the splittings 
obtained here should be viewed as a lower bound, since 
only effects of electronic disorder are taken into account, 
i.e., we assume that in the absence of such disorder the 
two states of the TLS are degenerate [^ . 

Using the distribution of splittings, we derive the tem- 
perature dependence of the zero-bias conductance of a 
metallic sample containing an ensemble of TLS. On this 
basis, we re-analyze the two-channel Kondo interpreta- 
tion of the zero-bias anomalies ||l7H19|] . The discrepancy 
between the observed scaling behavior and that derived 
in the present work presents, in our view, a significant 
challenge for the two-channel Kondo interpretation. Fur- 
thermore, the estimate of the total number of degener- 
ate (in the absence of electronic disorder) TLS which 
would be needed to produce the conductance observed 
in Ref. [|l7| is found to be unphysically large, indicating 
another problem with the two-channel Kondo scenario. 

The paper is organized as follows. In the next Sec- 
tion we consider in greater detail the role of symmetry- 
breaking disorder, present the main results, and discuss 
their implications for the interpretation of the Ralph- 
Buhrman experiments [[l7| . Section III contains the 
derivation of the splitting distribution for the magnetic 
and non-magnetic Kondo effects. A brief discussion and 



the summary are presented in Section IV 



II. DISORDER IN THE SPIN SECTOR. 

The general n-channel anisotropic Kondo Hamiltonian 
for a dynamic impurity f located at r = has the form 



a— 1 a 
3 



(1) 



where e is the Hamiltonian of band electrons, a^ (r) = 

ELi Ea/3 ^L (r) '^ifi'fPpa (r) is the electron spin-density 
operator at r, cTq,^ is the vector of Pauli matrices, and 
Jj are the exchange coupling constants. In the isotropic 
case we use the notation J = Jj. Greek indices are used 
to label spin quantum numbers, and Latin indices denote 
channel quantum numbers 




When considering the or 



FIG. 1. Schematic representation of atomic two-level sys- 
tem (TLS). 



convenient for our purposes to use, as a simple model of 
a TLS, an atom which can tunnel between two minima of 
a double- well potential located at r = ±b (Fig. |l|). The 
corresponding Hamiltonian has the form 

H = H, + ViT\ [4,^ (b) ^ (-b) + 4,^ (-b) ^ (b)) (2) 

+ ^3^3 (^t (b) 4, (b) - V^ (-b) i, (-b)) , 

where He describes free electrons, and ti and 1% are Pauli 
matrices operating in the two-state Hilbert space of the 
TLS. Vi is the amplitude of the "pseudospin-flip" process 
whereby the TLS tunnels between its two states accom- 
panied by the transfer of an electron from b to — b or vice 
versa. V3 is the position-dependent interaction between 
the TLS and electrons. The summation over physical 
spin is implied in the terms bilinear in electronic anni- 
hilation and creation operators ip ^nd t/"^. Terms pro- 
portional to r2 are absent in the Hamiltonian because of 
the combined effect of the invariance under time rever- 
sal and locality [|l|. Note that the apparently non-local 
term proportional to fi in Eq. ^ is an artifact of the ap- 
proximation neglecting the full momentum dependence of 
the coupling Vi. 

The Hamiltonian in Eq. (0) can be cast in the form 
of Eq. (|^) with corresponding couplings Q J3 = 
(27rpi?6/\/3) V3 and Ji = 2'kVi. The channel quantum 
numbers in Eq. dl^) would then refer to physical spin, 
while spin quantum numbers correspond to the impurity 
atom and electronic excitations being located at either of 
the TLS potential minima, or, depending on the choice 
of the basis, to different parity eigenstates. Anisotropy 
of the couplings is unlikely to occur in (hypothetic) mag- 
netic realizations of this model, while TLS realizations 
generically possess strong anisotropy |gl| . 

The disordered environment is modeled by adding to 
the Hamiltonian in Eq. (|l|) the term 



Dital two-channel realization 



of the Kondo effect due to electron-TLS interaction, it is 



iJdis =H, + H,= fdr [Uc (r) n (r) + U, (r) • & (r)] , 

(3) 

where n and & arc the charge- and spin-density opera- 
tors of the conduction electrons, and Uc and Us are the 
corresponding random potentials. Formally, the disorder 
Hamiltonian ifdis introduces two additional energy scales 
into the problem - the inverse scattering time t~^ and 
the inverse spin-scattering time t~^. 

Before proceeding further, we will comment briefly on 
the role of the charge disorder term He in Tfdis- It has 
been discussed in numerous works including Refs. [p|-p^. 
A related self-consistent model for the case of a finite 
concentration of Kondo impurities has been considered 
in Ref. |12|. This term affects only the charge degree of 
freedom, and one of its effects is to randomize the Kondo 
temperature Tk- In systems of lower dimensionality, it 
can also produce singular corrections to the energy de- 
pendence of thermodynamic and transport coefficients 
in the perturbative high-temperature regime. Incorpo- 
rating charge disorder into the description of the low- 
temperature (T < Tk) regime of the Kondo effect has 
not been achieved so far. However, it has been argued in 
Ref. gl that, in the one-channel magnetic case, the basic 
nature of the ground state as a local Fermi liquid would 
not change. 

In the context of the magnetic Kondo effect, we will 
only consider the role of spin disorder. Furthermore, we 
will restrict the consideration to the case of the ordi- 
nary one-channel effect, both because multichannel mag- 
netic realizations have not been unambiguously observed, 
and because the consideration of the renormalization of 



splittings (Section III) cannot be transferred to the case 
n > 2. [The renormalization of splittings in the n = 2 
case is discussed in Section UA. The results for a hypo- 
thetical magnetic isotropic n = 2 case can be obtained 
by a simple change 7m -^ 27^ in Eq. (g).] 

In the orbital two-channel realization of the Kondo ef- 
fect, the roles of spin and orbital degrees of freedom are 
partially interchanged. The potential scattering term He 
in Eq. (||) will produce a contribution analogous to the 
spin-scattering term, Hs, when Eq. (||) is transformed 
into the form of Eq. (|l|) [p3. In what follows, the term 
"spin scattering" is understood to apply either to phys- 
ical spin in the context of magnetic Kondo effect, or to 
"pseudospin" constructed out of electronic states of dif- 
ferent parity in the context of the orbital two-channel 
realization. 

Our choice of model for disorder assumes that break- 
ing of the spin symmetry is due to a well-defined set of 
scatterers present in the system. The scatterers couple 
to the spin degrees of freedom, e.g. "frozen" magnetic 
impurities in the magnetic one-channel case, or ordinary 
non-magnetic defects in the orbital two-channel case. Re- 
placing such a set of defects by a continuous random 
Gaussian-distributed potential, as is frequently done in 
transport calculations, is not warranted here because the 



distribution of splittings A is non-universal. That is, its 
form depends on the choice of the distribution function 
for the random potentials Uc (r) and Us (r) , and there- 
fore a more realistic model is required. 

In the model of isolated scatterers spin disorder is given 
by 



U,(r)=g^S,5(r-rO, 



(4) 



where Sj are randomly oriented frozen spins located at 
randomly selected points r^, and g is the corresponding 
exchange coupling constant. The distribution of each of 
Si is assumed isotropic and is given by 



VsiS) 



1 

2^ 



S{1 



s') 



(5) 



In the TLS case, we use a slightly more general expres- 
sion, allowing for scatterers of finite size: 



U, 



(r)=^[/(r-rO 



(6) 



Since "spin" degrees of freedom in this case are a subset of 
orbital degrees of freedom, the above expression contains 
both charge disorder and "spin" disorder terms. We will 
only concentrate on the effect of the latter. 

Below we will restrict our consideration to the case 
of rotationally invariant potentials [/, so that the cor- 
responding scattering matrix T can be reduced to the 
diagonal form T — diag(7;), where 71 can be expressed 
in terms of phase shifts r]i for each value of the orbital 
angular momentum as 



T, 



1 

■nv 



e**" sin?7;. 



(7) 



The coordinates r^ of TV impurities are drawn from a 
uniform distribution Vr ({ri}) — 1/V^ where V is the 
total volume of the sample. The calculations are per- 
formed in the limit N^V ^ 00 with the concentration 
c = N/V kept finite. It is assumed that the concentra- 
tion is small in the sense that the typical inter-impurity 
distance d ~ c~^'^ is much larger than the Fermi wave- 
length \p — 211 Ipp, which implies, in the magnetic case, 
EpTs ^ 1, and EpT :^ 1 in the TLS case, where Ep 
and pf are the Fermi energy and Fermi momentum, re- 
spectively. The effective "spin" scattering time in the 
TLS case is also of the order of r, and we will keep the 
notation Ts for it. 



A. Splitting of Internal States of Dynamic Impurity 
by Spin Scattering. 

What is the main effect of spin scattering on the be- 
havior of a Kondo impurity? Even though a non-zero 
value of T~ ""^ means that "spin memory" of electrons has 
a finite lifetime, the spin-flip processes at the location of 



the Kondo impurity still lead to logarithmic divergences 
in the high-temperature perturbative expansion (see Ap- 
pendix H). In fact, T~^ does not directly compete with 
the Kondo temperature Tk- Nevertheless, spin scatter- 
ing can change the low-temperature behavior of an im- 
purity by introducing an effective splitting. A, between 
its internal states. 

To understand qualitatively why finite t^^ does not 
by itself destroy the Kondo effect, consider the under- 
lying Anderson model for a magnetic impurity. In this 
description, the Kondo effect is reflected in the logarith- 
mic divergence of the perturbative contributions to the 
impurity electron self-energy at Ep. There are two pro- 
cesses that contribute to this self-energy: tunneling of 
an impurity electron into the conduction band, and the 
reverse process, tunneling of a conduction electron onto 
the impurity. Both of these contributions are logarith- 
mically divergent at Ep but with opposite signs. For 
an impurity without on-site interactions, the two terms 
cancel and there is no Kondo effect. Interactions remove 
this cancelation, essentially because an occupied site can 
only decay in one way - by an electron tunneling out 
(double occupancy is forbidden or strongly suppressed 
by Coulomb repulsion), while an unoccupied site can de- 
cay in two ways - by an electron of either spin tunneling 
in. The appearance of a finite t~^ may change slightly 
the relative tunneling rates for spin-up and spin-down 
electrons, but it cannot significantly change the factor of 
two difference between the rate of decay of an occupied 
and an unoccupied site. Hence t~^ does not directly de- 
stroy the logarithmic divergences in perturbation theory 
associated with the Kondo effect. 
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FIG. 2. (a) The dominant contribution to the effective 
splitting A. The solid line represents the conduction elec- 
tron Green function, and the broken line represents pseud- 
ofermions. The dot represents the interaction vertex J, and 
the cross corresponds to impurity scattering, (b) The lead- 
ing logarithmic order (parquet) contribution to the splitting. 
(c) The sum of parquet diagrams. The large dot represents 
the renormalized (in parquet approximation) scattering am- 
phtude Jr. (d) Renormalization of splitting. 



Importantly, however, the main effect of spin scatter- 



ing is to break time-reversal invariance and hence induce 
a splitting A between the two spin states of the Kondo 
impurity. The splitting results from the appearance of a 
random non-zero quantum-mechanical expectation value 
of the local spin density of states cr (e) at the impurity 
site. The coupling of cr (e) to the dynamic impurity via J 
is responsible for breaking the energy degeneracy between 
different orientations of the impurity spin. In the dia- 
grammatic language, the splitting is associated not with 
the modification of the standard set of logarithmic dia- 
grams, but rather with proliferation of a new set of dia- 
grams which were forbidden by SU(2) and time-reversal 
symmetries in the absence of spin disorder. The leading 
contribution to A in perturbation theory in the coupling 
strength J is shown in Fig. gja. The energy scale estab- 
lished by A serves as a cut-off of all logarithmic diver- 
gences in the perturbation theory in J. 

Formally, the main effect of spin disorder is to gen- 
erate a self-energy term which is essentially energy- 
independent, and can be reinterpreted as a contribution 
to the effective impurity Hamiltonian. It is a Hermitian 
2x2 matrix in the space of impurity states. In the mag- 
netic case, it is convenient to expand this matrix in the 
basis of Pauli matrices as 



H. 



eS 






(8) 



where the Aq, denote the components of the impurity 
energy in this basis. 

The discussion presented so far does not distinguish 
between the one- and two-channel cases. This is natural 
since the high-temperature diagrammatic expansions in 
the two cases have identical structures, and differ only 
in factors of 2 (from the two channels) for closed elec- 
tronic loops. If the Kondo temperature Tk < A, spin-flip 
processes are suppressed at T < A, the strongly corre- 
lated state is never formed, and below A all tempera- 
ture dependences are of the Fermi liquid type. If, how- 
ever, Tk > A, the behavior of the one- and two-channel 
Kondo systems is very different. In the one-channel case 
a non-zero but small A is equivalent to a weak polarizing 
field acting on the Kondo impurity, resulting in finite but 
small changes in the values of impurity susceptibility and 
conductance. However, since these quantities depend on 
Tk, which is itself altered in a random way by disorder, 
no significant experimental consequences follow from a 
small splitting A < Tk- 

Conversely, in the two-channel case when A 7^ there 
appear two distinct low-temperature regimes. A non- 
Fermi-liquid regime survives in the interval between T* < 
T < Tk, where T* is a new characteristic temperature, 
T* = A'^/27rTK [|l|-||. A Fermi-liquid regime emerges 
below T* . Hence, the low-temperature properties of an 
ensemble of two-channel Kondo impurities will explicitly 
depend on the distribution of splittings A. 

If A > Tk, a non- Fermi- liquid regime does not exist, 
and at temperatures below A the splitting becomes the 



only relevant energy scale. In other words, the TLS (or, 
rather, the composite object comprised of the TLS and 
correlated electrons) becomes frozen in its lowest energy 
state. Excitations above this state result in T^ depen- 
dences which blend with other Fermi-liquid effects that 
are always present. As was demonstrated by Moustakas 
and Fisher |Q , a generic set of TLS parameters almost 
always corresponds to A > Tr-. In order to observe the 
non-Fermi-liquid square-root scaling behavior, the pa- 
rameters of the Hamiltonian have to be fine tuned |3J] . 
The distribution of splittings in the ordinary (one- 
channel) magnetic case has a simple form (Ref . f|l6| , Sec- 
tion flT 



7'™(A) 



47„ 



(A2+7,2„ 



(9) 



where 7m = {"^.n^c/Spp) Jgu'^Ep is a constant which de- 
termines the scale for the typical values of A. It is pro- 
portional to the strength of the dimensionless Kondo cou- 
pling vJ , where v is the density of states at the Fermi 
level, and to the magnitude of spin scattering g [Eq. (Q)]. 
The quadratic suppression of Vm at small A results from 
the fact that A^ is a sum of three random variables A^ , 
each possessing a smooth distribution near A^, = 0. Cru- 
cially, the Kondo temperature Tk ~ Sfg"^/*-''''-' is expo- 
nentially small, so even in systems with weak magnetic 
disorder (c/pp, vg <^ 1) 7m may be comparable to Tk- 

In the orbital two-channel Kondo case, we will show in 
Section |l| that V (A) is described by a more complicated 
analytical expression. 



V{A)= / - — p^ 



d9 K{e) A-/ (6) 



1K3 [A2 + 72(6i)]^/^ 

where K (0) = {cos^ 9/Ki + sin^ e/K3y\ 1 [9) 

27rc 

3^ 



(10) 



27rc K{e) 

J— ^F {6), and 



F{e) = l 



cos 



In |tan( 



1 



Icosei 



(11) 



The constants Ki and K^, depend on the strength of elec- 
tronic coupling to the TLS and static scatterers, 

K^=2{EFViut)\ Ks^2iEFVsiytf{2pFb)\ (12) 

where the scatterer strength t is expressed in terms of 
the scattering phase shifts rji introduced in Eq. (M): 



t = ^(-l)'(2? + l)sin(27/i). 



(13) 



The function 7 (9) is the analog of the parameter 7m 
introduced in the magnetic case. 

In the asymptotic limit A 3> —r- max < WKi , WK^ > 

p%, L J 

the above expression simplifies to 



P(A) 



V2cK^K^ 
3p3,A2 



F {9) d9 



{Ki sin^ 9 + K^ cos2 9) 



3/2 ■ 



(14) 



The distribution function in Eq. ( p^ ) has the same 
asymptotic A"^ behavior as Eq. ^ for large A. How- 
ever, the full distribution function given by Eq. ( p^ ) 
is linear rather than quadratic at small A because, in 
the absence of T2 terms in the bare Hamiltonian, A^ is 
a sum of two rather than three random variables. In 
physical models of TLS the couplings are usually re- 
lated via Vx - Vs{2pFbf (see Ref. |2l[). This typi- 
cally corresponds to K^ ^ Ki, so that the distribution 
function acquires a third asymptotic region. It extends 
between the sharp maximum at A '^ (c/p^) \/Ki ~ 
(^vVict/p^^ Ep and the beginning of the A~2 decay at 
A ~ {c/Pf) VKz ~ {yV^ct/p^p^ Ef, and corresponds to 
approximately linear decay of the distribution function 
from a constant value. If Ki w K^, the intermediate 
asymptotic regime disappears. 

The graph of T-" (A) for the strongly anisotropic regime, 
corresponding to the choice of parameters discussed in 
Appendix g (c/p|, « IQ-* K3 = .38EJ,, Kx = 3.7 • 
IQ^^E^), is shown in Fig. g[ The isotropic case {Ki = 
K3 = .5Ep) is shown in the top inset. Note that both 
crossovers between asymptotic regimes in the main graph 
occur at values of A which are significantly smaller than 
the r.m.s. sphtting quoted in Ref. ||l|] (see also Ref. |p5|): 

(A2)^/^ - lOOK ~ IO-^Ef, where the value Ep ^ 8 ■ 
lO'^K for copper has been used. 




FIG. 3. The distribution function of splittings A. The pa- 
rameters of the distribution function are chosen as folIows(see 
Appendix^): c/p% = lO"*, K3 = .38EJr, Ki = 3.7 ■ 10'^ EJ,; 
the Fermi energy of Cu is Ef = 8 ■ lO^K. Insets: (a) Isotropic 
case Ki = K^ = .dEp. (b) Distribution of the crossover tem- 
perature r*; parameter regime is the same as in the main 
graph; Tk ^ 8.2K. 



Assuming that fluctuations of the Kondo temperature 
Tk can be neglected (see Appendix |c|), the distribution of 
T* - the temperature at which the crossover from non- 



Fermi-liquid to Fermi-liquid behavior occurs 
inferred from Eq. (|lC|): 



ViT*) = 



d6 K{9) -f{9)/VWr, 



K 



A-K ^KiKy, [T* + 72 (5)) I2t,Tk 



is easily 



(15) 



3/2' 



This expression is applicable only for T* <C Tk, where 
the relationship T* = ^^ JI-kTk holds. The graph of 
V (T*) is shown in the bottom inset in Fig. pi. The cor- 
responding limiting behaviors are 



27rc 



Tr/2 



dB 



F2 (B) ' 



7'(T*)cx < 



T* < (c/p|,)' min (Xi, X3) /Tk, AGjv (T) « ^CT^/^T (0) / V (xT) dx, (17b) 



V^^T^T^ Jo "" VKlK's 

(c/p|,)' max (Xi, i^a) /T^ <T* <Tk, 

(16) 



where x = T*/T, and T' (T*) is a normahzed distribu- 
tion. The function T (a;) reflects "freezing out" of some 
impurity degree of freedom, and is likely to decay expo- 
nentially at large x. Thus the convergence of the above 
integral at large x is provided by either T or T' depending 
on the temperature. 

It has been conjectured in Ref. M] that the observed 
T^/2 scaling of Gn may be consistent with the existence 
of a distribution of splittings because of an auto-selection 
process, in which only impurities with sufficiently small A 
contribute to conductan ce. S uch a scenario would imply 
that the integral in Eq. (17a) could be approximated as 



'^/^ A-3/2 (61) i^ (61) 
da 



In the intermediate regime (c/p 

v2 



Ki/Tk < T* < 



c/pyj" K^/Tk the distribution function V {T* ) is pro- 



portional to 1/\/T*. Both the intermediate and the far 
asymptotic regimes may not exist if the corresponding 
boundaries become comparable to, or larger than, Tr-. 



B. Relevance of Disorder to 

Two-Channel-Kondo-Model Interpretation of 

Zero-Bias Anomalies. 

If the parameters of TLS are randomly drawn from a 
certain distribution, the net contribution to conductance 
AGn (T) from N TLS can be approximated as 



AG 



N 



(T) ^N{G)+0 (\/7v) 



where (G) is the contribution of a single TLS averaged 
over the distribution of T* . It was estimated in Ref. |lj] 
that at least 10 separate TLS in the vicinity of the 
nanoconstriction, and up to 40 in some samples, must 
be present in order to explain the observed magnitude of 
the zero-bias anomaly. In this regime the second term on 
the right-hand side in the above expression would mani- 
fest itself as noise on the experimental curves, and hence 
can be neglected. A scaling ansatz for the anomaly in the 
zero-bias conductance due to a single TLS can be repre- 
sented aiT <Tk as AG(T) = CT^/'^T{T*/T), where C 
is a constant and T is a smooth function with the lim- 
iting behavior T(0) = 1 and T(a; ^ 1) -^ 0. The signal 
from N TLS is then written as 

AGAr(T) =iVGri/2 I 'p{T*)T[T*/T)dT* 

= NGT^^^ f V {xT) T (x) dx, (17a) 



resulting in AGa? = NGT^^"^ by virtue of the normal- 
ization of V. However, this approximation is only valid 
i f V (T*) is peaked so sharply that the integral in Eq. 
(17a) is not much different from its value in the limiting 



case V (T*) c>c S (T*). We will now consider under what 
circumstances such a behavior of V (T*) is possible, and 
what consequences a different behavior of V (T*) would 
have. 

Randomly formed TLS in metallic glasses typically 
have a broad distribution of as ym metries A^ even in the 
absence of electronic disorder pq] . Neglecting contribu- 
tions from random A^^ at first, we can write V (A) — 
l/Wz, where Wz is independent of A in the region of in- 
terest. The corresponding distribution of T* is V (T*) = 
y/nTK/2T*/Wz, and consequently AGat must display a 
linear in temperature behavior, 



AGn = NC, 



a^r$^<^'- 



Non-zero values of A^^ distributed with a width Wx < Wz 
would only exacerbate the discrepancy with the experi- 
mentally observed AGat ~ T^/^ behavior: specifically, 
the distribution of A would acquire a linear dip at 
A < Wx, yielding a flat distribution of T*, leadi ng i n 
turn to AGat (T) -■ T^/^ from the integral in Eq. dlTa] ). 
Thus in order for the Kondo effect in its orbital two- 
channel realization to be the cause of the observed T^/^ 
scaling, one must assume a set of nearly degenerate TLS, 
at least before disorder is taken into account. Classi- 
ness as a source of TLS in Ralph-Buhrman samples must, 
therefore, be ruled out. 

Let us now turn to the case when electronic disorder 
is the only source of TLS splitting, i.e. the TLS are 
assumed to be formed by some mechanism which, in the 
absence of coupling to conduction electrons, ensures their 
degeneracy. The analytic form for the distribution of T* , 
Eq. (15), can be substituted in Eq. ( |17a| ) together with 
an exponential ansatz T (x) = e~^ . Using the available 
experimental data to determine the parameters of the 
distribution V (T*) in Eq. ( [l5| ) is not straightforward. 



and is discussed in detail in Appendix g. Integrating 
over X we obtain the following integral representation for 
the TLS contribution to conductance: 



C rift 

AGn (T) = iVCTi/2 / Jl 
Jo 27r 



d0 K{9) 



1- 



7(g) 
V27V7 



exp 



7^^) 



2ttTkT 



erfc 



7(0) 



V2i7VT 



(18) 



The remaining integral over 9 is performed numerically, 
and the resulting graph for the temperature dependence 
of AGn is shown in Fig. H. As discussed in Appendix H, 



a 

< 




-9.0 



-4.0 



-3.0 



-1.0 



0.0 



-2.0 

ln(T/T,) 

FIG. 4. Temperature dependence of the change in 

conductance due to scattering by an ensemble of TLS, 
based on Eq. (|l8|). The solid line is the result of nu- 
merical evaluation, and the dot-dash line is the linear fit 
In {AGn/NC) = -4.83 + 0.84 In [T/Tk]- 



electronic disorder is assumed to be caused by strongly 
scattering non-equilibrium vacancies. The plot in Fig. |j 
does indicate a behavior close to a power law. However, 
the best fit in the region between .OSTjc < T < .^Tk, 
corresponding to the interval AK < T < AK for Tk ~ 
SisT, gives AGn oc T"«* with afit « .84. The deviation 
from the T^/^ behavior expected for degenerate TLS is 
rather significant. While the measured exponent in the 
temperature dependence of zero-bias conductance does 
show deviations from a ~ 1/2, they do not exceed .25 
(ReL Pol), so that the value .84 clearly lies outside the 
experimental error. 

This result can be understood better by examining un- 
der what conditions Eq. (17b) can be valid. Let us as- 
sume, for simplicity, that V is characterized by a sin- 
gle parameter, its width W . Formally, Eq. (17b) can 
be used when T ^ W . Let us consider a hypothetical 
case of a sufficiently sharp distribution, e.g. a Gaussian, 
V 



^2/ttW^ exp [-T*V2VF2]. Substituting this form 
into Eq. (17a) we find that in the temperature region 



W < T < lOVF, the best fit gives a power-law expo- 
nent of approximately .7, while the asymptotic behav- 
ior .5 is approached within a 10% accuracy in the re- 
gion IQW < T < 2QW. Thus, even in the case of a 
sharp distribution, the relation T <^W has to be under- 
stood as implying at least an order-of-magnitude differ- 
ence. Returning to the actual distribution Eq. dlQ), we 
note that it is a much broader function with a power- 
law 1 /r*3/2 ^g^j]^ gQ thdX the condition for the validity of 
Eq. (17b) is even stricter. At the same time, the larger 
of the two parameters controlling the width of V (T*) is 

W - {c/pI^Y Ks/Tk - 3K - .3Tk, so that the asymp- 
totic condition W ^ T can never be satisfied for tem- 
peratures below Tk- 

It should be observed that the most optimistic choice 
of parameters for the two-channel Kondo interpretation, 
corresponding to scattering by non-equilibrium vacan- 
cies, leads to an estimate for the concentration of de- 
fects c/p^p « 10""*. Using such a choice of parameters to 
support the two-channel Kondo interpretation also meets 
with a difficulty concerning the fraction of defects which 
form two-level systems. Indeed, using the estimate made 
in Ref . |lj of up to 40 separate TLS in the vicinity of the 
nanoconstriction, the correspondin g es timate for the den- 
sity of active TLS is found in Ref. ^^ to be 10~'*/atom, 



or Cactivo/p;? 



2-10 "".It can be assumed that TLS with 



T* < To, where To is the appropriately chosen cut-off 
temperature, are active. Then the ratio of the concen- 
tration of active TLS to their total concentration cxls is 
given by 



^active 
CTLS 



V{T*)dT*. 



(19) 



Choosing Tq ~ .8ii', which is a slightly generous assump- 
tion, since the T^'^ behavior is traced experimentally 
down to temperatures T w AK, we find Cactive/cxLS ~ 
0.4. The total density of TLS is consequently estimated 
at ctls/pI ~ (1/0.4) Cactive/PF « 0.5 • 10-"*. Compar- 
ing this to the above estimate of c/p\ « 10"''' for the 
total concentration of defects forces the improbable con- 
clusion that half of all the defects in the constriction are 
two-level systems. In other words, in order for the auto- 
selection mechanism to work, the total number of TLS 
which would be degenerate in the absence of disorder 
must be so large as to be inconsistent with the results of 
measurements indicating a rather small overall density of 
defects in the nanoconstriction. Although indirect, this 
reasoning serves, in our view, as another indicator of in- 
ternal consistency problems with the two-channel Kondo 
interpretation of zero-bias anomalies in Cu nanoconstric- 
tions. 



III. THE DISTRIBUTION OF SPLITTINGS. 



A. Local Spin in a Spin Glass Environment 

The coupling of the conduction electrons to the random 
spins results in the appearance of a non-zero expectation 
value of the conduction-electron spin density at a generic 
point in the sample. When this spin density is coupled 
to the dynamic impurity, the lowest order effect is to 
generate a self-energy matrix, whose components A^ in 
the basis of Pauli matrices and to the lowest order in 
impurity concentration c are (see diagram in Fig. 0a) 



^J9Y.j^G\e■v.,)Sf, 



(20) 



where the zero-temperature Green function in the coor- 
dinate representation is given by 

G{e;v) = -^!^expii {pFr+—r] Sgne I . (21) 
ppr y \ VF J J 

Note that A^ is proportional to the RKKY-induced ran- 
dom spin polarization in direction a at the position of 
the impurity. To the lowest order in J, the distribution 
of splittings follows the distribution of internal random 
magnetic fields in spin glasses. The latter has been de- 
rived in Ref. [ [L6[ using a somewhat simplified RKKY in- 
teraction, in which its oscillatory character is modeled by 
random signs. The derivation below, while reproducing 
the essential results of Ref. [Q serves primarily to in- 
troduce the more technically involved derivation for the 
non-magnetic two-channel case presented in the next sub- 
section. 

Integrating over energy we obtain 



A„ = Jgv^ V ^^^ cos {2pFn) S^. 

J KPFn) 



(22) 



The distribution function of A = ^J'Yla ^"a therefore 
takes the form 



N 



N 



n.(A,=2A/nt/nt 

i i 

X <5 (1 - S,^) W A^ - :| ^ / (r,) / (r,) S,;S, J , (23) 

where K = 2 {aJg^-'Epf and / (r,) ee /, = ^"'^^^^3'^ 

ippri) 
Exponentiating the second (5-function and introducing a 
shorthand notation F = ^. /^S^ we rewrite the distribu- 
tion function as 

n„(A)=2A/n^/nf^(i-sf) 



2^' 



.ifil^^-^tiKF^ 



(24) 



Decoupling the last term in the exponent with the help 
of a "Hubbard-Stratonovich" transformation we obtain 



N 



V...{A)=2Ajl[f^S{l-S^)J-J^e^^^-' 



dfi 
2^' 



dX 



= 2Al%^'^^' 

27r 



f zA2 

d\ 



iXF 



[2m^K) 



3/2 



''''P I 2^x1 



/|f.a-s») 



-i/(r)AS 



N 



(25) 



The expression in the square brackets can be transformed 
as follows: 



/y/f^(l-S^)exp{-z/(r)AS} 



N 



drsin[|A|/(r) 



-,N 



exp 



V |A|/(r) 



dr 1- 



sin[|A|/(r)] 
|A|/(r) 



(26) 



To compute the last integral the following approximation 
is employed: since the dominant contribution to the dis- 
tribution function is expected to come from r ~ d, it is 
possible to decouple fast oscillations in / (proportional to 



cos 2pFr) from the slow decay {ppf) ■ Formally, / (r) 

1 1 • 1 /■ / \ cos f . , . , 
IS replaced with / [if, r) — ^ with a simultaneous 

(ppr) 

[dip 
dv I — . Integration over r now 

27r 



replacement I dv - 
gives 

dip TT^C |A| |cos(^ 



exp 



27r 



ip% 



exp 



27rc|A| 
■ipl 



(27) 



Absorbing ^J\iji\K into A we can rewrite the distribu- 
tion function as 



V^iA) = 2Aj^e^'^^' 



An 



(27riSgn/i) 



3/2 



A^dAexp ■ 



iSgn/i 2 2TTCy^\p\KX 



3pl 



(28) 



Rotating the contour of integration by (7r/4) Sgn/i we ar- 
rive at 



^m(A) 



^g»MA= 



2tt V tt 

27rc 1 + iSgn/i 
" 3^ V2 



2 /■ 2 , f 1 2 

X dxexp ■^ — X 



VMkx 



(29) 



Performing another change of variables y^ = 2 |/x| A^ 
and rotating the contour of integration over y we obtain 



the following integral representation for the distribution 
function: 



X expl~-{x'^ +y^) +i^xy\ , (30) where h^ ^ h{r. 






Kifjj + Ksh^h, 



(35) 



27rc /ii: 



where 7™ = 7r^\ ~:r- Using polar coordinates {p,0) in 
opp V 2 

the (x, y) plane and integrating over p the distribution 

function can be rewritten as 



^ ,,, 2^ d r'^ cosBdO 

Vr,, (A) = Re, 



TT d^^J, ^,_^^^i^2^^) 



3/2' 



(31) 



The remaining integral can be performed by elementary 
means leading to the expression for the distribution func- 
tion (cf. Ref. PI) quoted in Eq. (|). 



B. Non-Magnetic TLS in the Presence of Static 
Defects 

The electronic contribution to the splitting between 
the energy levels of TLS is determined by the difference 
between the two eigenvalues of the self-energy matrix fll] 



Vi5-+ -V^5- 



(32) 



The components of this matrix are expressed in terms of 
the scattering matrix T corresponding to the potential 
U (r) introduced in Eq. (ph as 

5++ = -2iY,j ^jdvjdv'G,{v + h) 

X r, (r - r„ r' - r,) G, (r' + b) , (33a) 

and similarly for 5 The off-diagonal term 5-\ is 

6+_ ^ -2iY, j ^ j dv j dv'G,{v + h) 

i 

x%{v-v,y -v,)G,{v' -h). (33b) 

The factor of 2 in front is due to summation over channels 
- the two orientations of the real spin of the conduction 
electrons. The resulting electronic contribution to the 
splitting is 

A^ = Vi [Re {5++ 5^^)f + AV^ (Re<S+_)^ . (34) 

After integrating over e and approximating jr^ ± b| = 
r,; ± b • rii, rii = Vi/vi, the distribution function for A is 
given by an expression analogous to Eq. (p2|): 



sin2pi^rj . .„ , . , 
g- sin(2pi?b -Hi), and con- 

stants Ki, K3, and t were defined above in Eqs. (y_2|,^_3|). 
K3 is related to K3 via K3 = {2pph) K3. 

Following the technique used in the preceding section, 
the distribution function can be represented as an in- 
tegral over a Lagrange multiplier /i and two "Hubbard- 
Stratonovich" variables Ai and A3, 



V{A)^2A I '^e'^^" 



dXidXs r Xi 

exp < I 



2T:ip\/ K1K3 



X2 
^^3 



2nK3 3pj, 



— 3-m(Ai,A3) 



2fiKi 
(36) 



where the function u is defined by an expression analo- 
gous to Eq. (|^), 

47rc f 

— 3-u(Ai,A3) = c / dr[l-cxp{-zAi/(r) -iXahir)}]. 

(37) 

Integrating over orientations of n, and decoupling fast 
oscillations in / and h, we arrive, in the approximation 
sin2pi?6 « 2pi?6, at an analogue of Eq. (p7|): 



u(Ai,A3) 



dC f^'' dip 



e .In 2^ 



1- 



sin I A36sin<y9C J 



X3bsm(p(^ 
X exp{— lAi cos(/jC} 7 (38) 



where b = 2ppb and C = 1/ ijapr) . Integration over ip 
now gives 



m(Ai,A3) 



26 I A3 



=dC 



|Ai|sinCj^_^/^C|Ai| 



I A3 1 6 



6 I A3 



(39) 



where Ki^ is the modified Bessel function. Completing 
the remaining integration over C, and rescaling 6A3 — > A3 
so that K^ in Eq. ( |3^ ) is replaced with K^ = IPK3, we 
obtain 



u [x, y) = ^ In '^' + Yf ^^' + v/^^+^- (40) 

\y\ fI 

Using polar coordinates (A, 0) in the (Ai,A3) plane, in- 
troducing another set of polar coordinates (p, ■(/') in the 
(A, ji) plane and integrating over p we obtain the follow- 
ing integral representation for the distribution function: 
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Jq t^ 



"" dB K{9) d /""/^ d^p 



^Ay/K^ d-fJo ijsm2i;-l 



(41) 



The final result obtained after integrating over tp, to- 
gether with the definition of 7 (9) , has been presented in 
Eq. ®. 



C. Higher Order Contributions. 



The scaling ansatz used in Eq. (17a) depends crucially 
on the shape of the distribution function V in the re- 
gion of small splitting A. Thus, the applicability of the 
above analysis hinges on whether the perturbative cal- 
culation of A in the preceding subsection is sufficient 
for A's smaller than Tk- Only the lowest order diagram 
in Kondo coupling has been retained in the calculation 
so far, and we now turn to the consideration of higher 
order contributions. 

These higher order contributions can be conveniently 
separated into two classes. The diagrams belonging to 
the first class (Fig. Ep) correspond to the renormaliza- 
tion of the Kondo scattering amplitude with parquet dia- 
grams, which collect the leading order logarithmic terms. 
They produce corrections to A which are smaller by pow- 
ers of the Kondo coupling, and, importantly, unlike the 
Kondo scattering amplitude itself, do not contain any 
logarithmically divergent terms. These contributions can 
therefore be neglected. 

Indeed, the contributions of the diagrams belongingto 
the first class is typified by that of the diagram in Fig. gb: 

SA^ =. (lyg) {vjf ^ -^^^ cos {2ppn) ^f 

. {VFTi) 

In -^ + Const. (42) 

This expression is valid as long as 

vp/n < Tk, (43) 

which, for Tk/Ef ^ 10^** would be violated only at un- 
physically small concentrations c/p^ < 10~^^. It is seen 
that this contribution does not contain any uncontrolled 
logarithmic divergences. The reason is that the sum of 
all diagrams of this type can be written in the form of Eq. 
(P0|), where the bare scattering amplitude J is replaced 
with the renormalized amplitude Jjj(e, e). The renor- 
malized amplitude is singular at small energies, but the 
singularity is in teg rable, and after the integration over 
energies in Eq. (EG) it only manifests itself in finite loga- 
rithmic terms like the first term in square brackets in Eq. 
(P^). The inequality (43) ensures that the contribution of 
these terms is small in the parameter (i/ J) In ^^^^ ^ 1, 
and can be ignored. 



In contrast, the diagrams of the second class corre- 
spond to a further renormalization of the scattering am- 
plitude by particle-hole pairs, and retain the usual log- 
arithmically divergent factors. These terms were previ- 
ously analyzed perturbatively using the renormalization 
group (RG) approach in Ref. [^, where it was estab- 
lished that in anisotropic models the effect of these terms 
may be to renormalize downwards the effective splitting 
at the scale just above Tk- We argue that this anal- 
ysis cannot be straightforwardly extended through the 
crossover region into the low-temperature regime. In- 
stead, we show that the effective splitting at tempera- 
tures below the Kondo temperature can be deduced from 
the lowest order perturbative result on the basis of uni- 
versality properties of the two models under considera- 
tion, the one- and two-channel Kondo models. We em- 
ploy the language of the underlying Anderson model, as 
it affords a unified description of the high- and low-energy 
regimes [g[. 

The simplest diagram of the second class is shown in 
Fig. ||d. These graphs correspond to linear (in A) screen- 
ing of the splitting by the Kondo interaction. In the per- 
turbative renormalization group analysis by Vladar and 
Zawadowski p^ it was found that the splitting tends to 
be renormalized downwards, and in the case of strong 
anisotropy, at least one component of the splitting may 
be renormalized significantly. 



A,,eff(T>TK)- A, 



/ vJi 



\AlyJ3 



l/iuJs 



However, the perturbative analysis in Ref. |2l| cannot be 
extended to energy scales below Tk- An extrapolation of 
the perturbative RG results from the energy scale E > 
Tk to predict the values of A at T < Tr- (Ref. ||§) is 
unjustified because splitting is a relevant operator in the 
RG sense. 

Let us consider first the isotropic one-channel magnetic 
case. The perturbative renormalization group calculation 
of Ref. [|l| , although formulated in terms of the orbital 
two-channel Kondo model, can be transferred, with mi- 
nor modifications, to the one-channel case as well because 
of the identical structures of the corresponding high- 
temperature perturbation expansions. When graphs of 
the type shown in Fig. ||b are neglected, the splitting has 
the literal meaning of an external magnetic field h act- 
ing on the impurity. The unrenormalized impurity Green 
function has the form 



g ^ {lo - ed + ii-- (^y 



(44) 



where e^ is the energy of the singly occupied impurity 
state. The ground state of the model can be described as 
an effective Fermi liquid, in which the Green function of 
the impurity spin retains its form under strong renormal- 
ization. The fully renormalized Q will contain additional 
self-energy contributions [@|, 



Ll) — ed + iT/z + /i • cr — E [uj, h) 



(45) 
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where F is the width of the Kondo resonance. The self- 
energy E = Eq + S ■ (T is expanded at small lu and h 
as 



Eo {uj,h) = Eoo + 

s = e;>, 



1- - I w + iO(w^ 



(46) 



where EJ^ and z are the renormalization constants, and 
Eqo is the constant term which, in the case of symmetric 
Anderson model, is equal to ej, ensuring that the Kondo 
resonance is centered at the Fermi energy. 

Factoring out the quasiparticle weight z, the effective 
field heff takes the form 



h,ff = z{i-j:'jh. 



(47) 



The perturbative RG calculation of renormalized split- 
ting in Ref. pW is equivalent to computing the same 



prefactor R (uj) = z (1 — EjJ = 



1 - dY./dh 



at a finite 



1 - dH/du 

frequency oj (see Eqs. (3.5)-(3.10) in the second of 
Refs. [pil). The frequency dependence of both deriva- 
tives in this region is logarithmic, so that, with logarith- 
mic accuracy, frequency can be identified with energy 
scale in the RG sense. Such scale-dependent quantities 
determine the properties of the system at temperatures 
of the order of the energy scale. 

The perturbative RG analysis cannot be continued be- 
yond some intermediate scale Eq > T^, and its usefulness 
for determining low-temperature properties relies on the 
assumption that further renormalization from E '^ Eq 
to E ^ Tk does not change the value of the renormal- 
ized quantity in any essential way. This assumption may 
be violated when renormalization of relevant operators is 
considered, as is indeed the case in the models considered 
here. 

At zero temperature, the prefactor in Eq. ( p7| ) is just 
the universal Wilson ratio R (Ref. |g]). Substituting its 
known value, i? = 2, we obtain a seemingly counterin- 
tuitive result that, despite downwards renormalization of 
splitting aX E > Tx, the splitting is actually increased at 
E ^ Tk by a factor of 2 compared to its bare value. Of 
course, at zero temperature the role of the weak effective 
magnetic field acting on the impurity is to polarize the 
Kondo screened complex (as long as A <C Tk), and hcs 
induces a "splitting" only in this sense. 

To explain the non-monotonic behavior of the effective 
splitting as a function of energy scale (or temperature), 
one should note that, in the Anderson model, the renor- 
malization of the splitting is proportional to the ratio of 
the impurity magnetic susceptibility x to C/T, where 
C is the impurity specific heat. This ratio is a non- 
monotonic function of temperature near Tk, essentially 
because of the maximum in the temperature dependence 
of Cat T-Tif. 

The Anderson model is less well suited for discussion 
of anisotropic Kondo models. Nevertheless, anisotropy 



can be modeled, at the price of additional potential scat- 
tering in the corresponding s-d Hamiltonian, by intro- 
ducing spin-dependent renormalization constants Zq and 
E'^o. Coupling anisotropy is an irrelevant operator (in 
both one-channel and two-channel spin- 1/2 Kondo mod- 
els, see, e.g. Ref. [Q), and the invariance of the Wilson 
ratio demands that 

2 = ^T (1 - s'tt) = ^i (1 - s'a) = v^i^i (1 - ^i) ■ 

(48) 

The choice of only two different renormalization con- 
stants Zf and zi corresponds to Ji = J2 ^ J3 in the 
Kondo effective Hamiltonian. If the field now is chosen 
in the xy plane (h-n in the above notations), the corre- 
sponding value of the self-energy at zero temperature is 
still 

huos = h^iy^z^zi (1 - E'lJ = 2h^i. 

The temperature dependence of this self-energy may be 
quite non-monotonic, as the transverse susceptibility de- 
viates strongly from the free-spin value at temperatures 
above the specific heat maximum. 

Let us turn now to the two-channel Kondo case. The 
above analysis cannot be transferred verbatim because 
the ground state is not a Fermi liquid, and T* does not 
have the meaning of the renormalized Fermi-liquid quasi- 
particle self-energy ]39[ |. In particular, the renormaliza- 
tion factor R (uj) can no longer be identified with the 
Wilson ratio of an Anderson-like model. Nevertheless, 
the salient features of this analysis survive. Once again, 
we can identify the self-energy contribution of Fig. 0a 
with an external field h acting on the impurity. The 
effective field in the high-temperature regime is renor- 
malized downwards, in the anisotropic case strongly pi| . 
However, in the low-temperature regime the weight of 
the quasiparticle excitations for which this effective field 
represents the self-energy is zero, so that this self-energy 
term does not define any physical low-temperature en- 
ergy scale. It has been demonstrated recently Ig^jQ that 
impurity thermodynamics in the low-temperature regime 
can be described in terms of three vector and one scalar 
Majorana quasiparticles. An external field h generates a 
self-energy contribution for the scalar Majorana fermion 
Em = h^ /2t:Tk = T* which is universal apart from its 
dependence on Tk fS^-Q ■ The universality [|6| ensures 
that T* is controlled by the unrenormalized value of h, 
namely the bare splitting A. 



IV. DISCUSSION. 

The results of our model calculation of the split- 
ting between the states of magnetic impurities and non- 
magnetic two-level systems (TLS) in the presence of spin 
disorder can be summarized as follows. First, the main 
feature of the distribution functions of the splittings is 
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their strong asymmetry. Formally the distribution func- 
tions in Eqs. (§) and ([ll| ) do not even possess finite first 
moments. This is a manifestation of the fact that all the 
moments of the splittings are dominated by disorder con- 
figurations with one or more scatterers located very close 
to the dynamic impurity. In real systems the shortest 
possible separation between the dynamic impurity and 
the nearest scatterer is determined by the lattice spac- 
ing. Therefore one has to introduce a short distance, lat- 
tice cut-off into the coordinate integration in Eq. (26). 




In the presence of such a cutoff all the moments are 
dominated by distances of the order of the cut-off, and 
thus, for example, the average rms splitting y (A^) is 
of the order of Ep [Jgiy^) \/c/pp for the magnetic case 
or Ep (max{Vi, V^ppb}) tvy^c/p'^ for the non-magnetic 
TLS case. These values are larger than the typical ones 
by a factor of \/p'p/c '-^ [ppd) » 1. This confirms the 
result obtained independently by D. Cox |El|. 

Second, since the splitting is the difference between two 
eigenvalues of a random Hermitian 2x2 matrix, at small 
splittings one observes the equivalent of level repulsion 
leading to a vanishing probability density to observe zero 
splitting. In the non-magnetic TLS case the matrix is real 
and symmetric, and the suppression is linear rather than 
quadratic. 

The asymmetry of the distribution functions has im- 
portant implications for the study of the effects of spin 
disorder on the behavior of Kondo impurities in the 
crossover and low-temperature regimes. Because of the 
build-up of many-particle correlations at low tempera- 
tures, the effective splitting between the levels of mag- 
netic impurities or TLS, Aefj (T), can be reduced at 
intermediate temperatures T ~ Tk compared to the 
bare value A (Refs. [gl],|0|). Since the distribution of 
splittings is very asymmetric, with the root-mean-square 
splitting much larger than the typical value, a proper 
renormalization analysis would necessarily treat the full 
distribution rather than just the first few moments. To 
what extent the distributions found here preserve their 
shape under renormalization is an open question. 

Nevertheless, the zero-temperature behavior of the 
splitting, or, more precisely, of the corresponding self- 
energy terms, is dictated by the universality properties 
of the one- and two-channel Kondo models, and can be 
extracted directly. This, in turn, made it possible to 
derive a corresponding temperature dependence of con- 
ductance for a collection of two-level systems. We find 
(cf. Ref. |11|]) a T" behavior with a w .84 in contrast to 
the T^/^ dependence observed experimentally p7| , p^ , pO{ . 
The concentration of the TLS which must be degenerate 
(in the absence of disorder) in order to sustain the two- 
channel Kondo interpretation is also found to be unphys- 
ically large. Both these arguments suggest that the two- 
channel Kondo model does not provide a consistent inter- 
pretation of the zero-bias anomalies observed in Ref. [0 . 

In conclusion, we have computed the distribution 
functions of splittings of magnetic impurities and non- 



FIG. 5. The lowest order logarithmic contribution to sus- 
ceptibility in the magnetic Kondo problem. The wavy lines 
correspond to the external magnetic field. 



magnetic two-level systems (TLS) induced by disorder 
scattering of conduction electrons. In the magnetic case 
these splittings only appear if the disorder breaks time- 
reversal symmetry, i.e. if the disorder is itself magnetic. 
In the non-magnetic case the degeneracy between the lev- 
els of a TLS is due to geometric symmetry about its cen- 
ter, and therefore is strongly broken by any type of dis- 
order. We find that the probability distribution of split- 
tings vanishes as a power law at small splittings, mak- 
ing nearly degenerate impurities a rarity. The typical 
values of splittings are found to be smaller than the av- 
erage estimated previously in Ref. |ll|. However, even 
in quite clean systems such as the ones studied in the 
experiments of Ref. |1^ , the broad profile of the distri- 
bution of splittings results in a temperature dependence 
of the conductance which is substantially different from 
the square-root law expected in the absence of disorder. 
Consequently, experimental observation of the square- 
root temperature and voltage dependences may not be 
a reliable indicator of two-channel Kondo physics. 
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APPENDIX A: 

To illustrate the small effects of spin disorder on the 
standard logarithmic terms in the perturbation-theory 
expansion, we consider the well-known lowest-order log- 
arithmic contribution to the spin susceptibility given by 
the diagram shown in Fig. o. 

As in other similar logarithmic diagrams, there appears 
a set of terms due to the off-diagonal in spin index parts 
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of the Green function. The corresponding analytical ex- 
pression is 



5x{T)_.P 



Xo 



deidt2 



3p(ei)/j(e2)-^crj (ei)crj (£2) 



1 — tanh -^p tanh ^ tanh -^ — tanh -1^ 

21 2J_ cyrp 21 21 



(ei -£2) 



(ei - £2)^ 



(Al) 



where Uj (e) is the electron spin density of states at the 
impurity site, and xo is the susceptibility of a free spin. 
Since the average {(Jj) = 0, there is no positive definite 
contribution to the logarithmic integral similar to that 
coming from the first term 3p(ei)p(e2)- Moreover, ab- 
sence of diffusive behavior for aj (due to the absence 
of an equivalent of the particle-number conservation law 
which enforces the universal diffusion pole in the density 
correlations) leads to the absence of "fine structure" in 
the correlator {(JjUj) at scales smaller than t~^. In the 
model of isolated impurities the effect of the second term 
in Eq. ( Al) reduces to a small correction to the coefficient 
of the leading logarithm E^ . 

Indeed, in the lowest order in concentration of mag- 
netic defects Sq, the local spin density of states is given 

by 



-(£) = -^ 



/ TIV 



KPFTa 



Ime-2^(P^+^/''^)''''S„, (A2) 



where Ya are the coordinates of the defects. Using 
{S^SD = -SabSij, we find 



{<y^{el)a,ie2)) =^ [dr(^ sin [2 
TT^ J \PFrJ 



sin [2 {pF + e2/vF) 



IT V {t^Qv) 



2 _C_ 

Pf 



ei +e2 
Ef 



ei -£2 



Ef 



(pf + f^i/vF)r] 



(A3) 



Substituting the last expression into Eq. (Al) we see 
that the second and the third term in square brackets 
produce irrelevant constants and terms small as T/Ef, 
while the first term results simply in a correction to the 

47rc 4 2 

5p% 



effective density of states b(? Iv^ — — tt^t'''^ (5^) 



APPENDIX B: 



The Kondo temperature in the TLS model is given 

by[il 



TK^EF^/{vJ^)(yh)\^^ 



J I \ ""^3 



(Bl) 



where Ji and J3 are electron-TLS coupling constants, 
and V is the local conduction electron density of states. 
Assuming a two-channel Kondo interpretation of the 
anomalies, the Kondo temperature is estimated in 
Ref. ijl^ to be between 5K and lOK. The choice of the 
bare coupling values vJ^ ss .2 and vJi « .007, which 
is adopted in the calculations leading to the graphs in 
Figs. 1^ and^, corresponds to the Kondo temperature of 
Tk ~ 8.2K. The respective values of the couplings V3 
and Vi in Eq. (||) are V3 « .74, and Vi « .0012, corre- 
sponding to the dimensionless distance between the TLS 
minima set at 2pFb « .15. 

The nature of defects in quenched vacuum-deposited 
films is not well understood, and is likely to vary de- 
pending on the details of a particular experimental setup. 
As a rough measure of disorder, the transport mean free 
path near the opening of the constriction is estimated in 
Ref. p9| | to be ^tr ^ 30nm in unannealed samples. The 
mean free path is related to the scattering phase shifts rji 
and concentration of defects c via Ga 



PF^r = 



Pf 



1 



4^"= Ei=i ' sin^ ('^i -m-iY 



(B2) 



where I denotes angular momentum channels. The trans- 
port mean free path in annealed samples is shown in 
Ref. ||l£] to be close to 300nm, i.e. it increases by a 
factor of 10 upon annealing. Therefore, most of the dis- 
order in the constriction is likely to be caused not by 
substitution impurities which cannot anneal, but by lo- 
calized structural inhomogeneities in stressed deposited 
films. Spatially extended defects, such as dislocations 
and grain boundaries, are unlikely to be a major cause 
of scattering for the same reason. 

The coefficients Ki and it's in Eq. (|l^ cannot be ex- 
pressed in terms of the transport mean free path alone. 
Their computation requires the knowledge of the phase 
shifts rji [which would ma ke it possible to infer the con- 
centration from Eq. (B2)]. Even allowing for an inde- 
pendent measurement of defect concentration (one not 
available for the Ralph-Buhrman samples), information 
about phase shifts is still needed if more than one of them 
is non-zero. 

It is seen from Eq. (|o|) that the scale for the typi- 
cal values of A is set by the combinations {c/pp) V^ii 
i = 1, 3, and, for small 77;, it is linear in both c and rji. To- 
gether with Eq. (B2), this implies that, for a fixed value 
of iti-, smaller values of 77; result in a broader distribution 
function V (A). For example, the choice of jsinTyj — 1/6 
leads to typical values of T* exceeding Tk- Conversely, if 
scatterers are strong (rji ~ 1), fewer of them are needed 
to produce the same value of £tr, and the distribution 
function is narrowed, implying smaller typical disorder- 
induced splittings. 

It should be noted that because of the Friedel sum rule, 
small values of all sin 77/ are possible only in the case of 
neutral defects. If defects are charged, at least one of the 
phase shifts cannot be small. 
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A typical neutral defect is a combination of a self- 
interstitial and a neighboring vacancy (a so-called Frenkel 
defect). We are not aware of any studies of whether such 
defects are more or less common in vacuum-deposited 
films than simple non-equilibrium vacancies, which are 
charged, and thus scatter strongly. However, the forma- 
tion energy of such neutral defects in Cu is 2 to 3 times 
larger than the formation energy of a vacancy ||44[| . Con- 
sequently, we concentrate on the non-equilibrium vacan- 
cies as the dominant source of scattering. Since one of the 
aims of this work is to explore whether the two-channel 
Kondo interpretation of the Ralph-Buhrman zero-bias 
anomalies can be compatible with the presence of the 
static scattering in their samples, the assumption that 
the bulk of disorder is caused by strongly scattering non- 
equilibrium vacancies is also quite appropriate. It leads 
to the smallest values of splittings, and is, therefore, the 
most favorable for the two-channel Kondo interpretation 
of the experiments in Ref . [|7| . 

Scattering by vacancies in Cu has been studied by sev- 
eral techniques (experimental, numerical, and their com- 
bination) in Ref. B^. The results of the fitting proce- 
dure combining experimental and numerical techniques 
quoted in Ref. [Esl give for the differences between the 
vacancy and the host phase shifts the following set of 
values: 770 ~ .92, 771 « —.7, and ry/>2 = 0. 
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dv{e) = i/(e) -z/= - — 






where r, are the coordinates of random scatterers. Sub- 



stituting the above expression into Eq. (C2) we obtain the 
following condition for the disorder-induced shift STk'- 



In- 



Tk + 6T, 



K 



Tk 



E 



tsa\{2pprj) , /2Tifr. 

{PFTjf 



VF 



0. (C4) 



The argument of the integral cosine ci is small as 

(Tk/Ef) [Pf/cY'^ -- 10"^, where the values Tk/Ef « 
10^"* and c/pp « 10~^ are used (see Section |IIB| ). 
It can therefore be approximated as ci (2TKrj /vf) ~ 
— C— In {2TKrj/vF), where C is the Euler constant. Typ- 
ical values of STk/Tk are thus seen to be of the order 

„3N-l/3l 



{Tk/Ef) [c/p],) 



and we estimate 



oit{c/pFf'^\n 

6Tk/Tk ^ 10^^. Therefore random variations of T/f due 
to static disorder can be safely neglected for the param- 
eter regime of Ref. [Il^ . 

Intrinsic variation of Tk among different randomly 
formed TLS is, of course, a possibility. However, it has 
been established in the main text that all TLS have to 
be degenerate in order for the two-channel Kondo inter- 
pretation to have a chance of succeeding. Thus one must 
assume that the TLS are formed by the same non-random 
mechanism, again excluding large variations of Tk- 



Fluctuations of the Kondo temperature Tk due to 
static disorder can be addressed in the same framework 
used for the calculation of the distribution of splittings. 
Indeed, in the leading logarithmic order, the Kondo tem- 
perature can be inferred from the scaling relation [pi| 



¥W0 



= -d\nD, 



(CI) 



where D is the renormalizcd band cut-off, i/) is a homoge- 
neous degree one function of dimcnsionlcss couplings J^i^, 
and <I> is a homogeneous function of degree two which de- 
pends on ip and bare values of the couplings. The effect 
of disorder is contained in the energy-dependent density 
of states V (D) . The homogeneity of both ip and $ al- 
lows one to make a substitution ip {Jv) —^ ip {J^)i where 
v is the average density of states, while transferring the 
energy dependence of v to the r.h.s. of Eq. (CI). The 
implicit equation for Tk takes the form 



d^ 

¥W0 



-Do 



Tk 



[v [D) + V (-£>)] dD 
WD ' 



(C2) 



where Dq ~ Ef is the unrenormalized bandwidth, and 
the boundary condition ip [J (Do)J_= was used. Follow- 



ing the consideration in Section HI B , the correction to 
the density of states as a function of energy e is expressed 
as 
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